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Abstract. We consider a multiferroic chain with a linear magnetoelectric coupling 
induced by the electrostatic screening at the ferroelectric/ferromagnet interface. We 
study theoretically the dynamic ferroelectric and magnetic response to external 
2 ' magnetic and electric fields by utilizing an approach based on coupled Landau- 

Khalatnikov and finite-temperature Landau-Lifshitz-Gilbert equations. Additionally, 
1-^ , we compare with Monte Carlo calculations. It is demonstrated that for material 

5IJ ' parameters corresponding to BaTiOs/Fe the polarization and the magnetization are 

O . controllable by external magnetic and electric fields respectively. 
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1. Introduction 

Magnetic nanostructures are intensively researched [1] due to their versatile use in 
technology. Multiferroics, i.e. systems that exhibit a coupled ferroelectric and 
ferromagnetic order, have received much attention recently [21 [31 H]. A variety 
of potential applications of multiferroics rely on a possible control of magnetism 
(electric polarization) with electric (magnetic) fields due to the magnetoelectric 
coupling. Indeed, recent experiments have demonstrated the existence of both 
effects[5l [6l [7]. Furthermore, a number of interesting phenomena associated with 
the coupled polarization/magnetization dynamics at interfaces and in bulk have 
been reported such as the electrically controlled exchange bias [8], the electrically 
controlled magnetocrystalline anisotropy[9] and the influence of electric field on the 
spin- dependent transport |10j. Theoretically, recent Monte-Carlo calculations for a 
three-dimensional spinel lattice were performed to study the magnetic-field induced 
polarization rotation [HI [12] . 

In this work we investigate theoretically and numerically the field-driven dynamics of 
the electric polarization and the magnetization of a ferroelectric/ferromagnetic system 
that shows a magnetoelectric coupling at the interface. For this purpose we consider a 
two-phase multiferroic chain consisting of 50 polarization sites and 50 localized magnetic 
moments, as sketched in FiglH The ferromagnetic (FM) part of the chain is a normal 
metal (e.g., Fe), whereas the ferroelectric (FE) part is BaTiOs (Fig. [1]). Recently 
this system has been shown to exhibit a magnetoelectric coupling [131 [H] and has 
been realized experiment ally [15] . The multiferroic coupling arises as a result of an 
accumulation of spin-polarized electrons or holes at the FE-insulator/FM-metal interface 
when the FE is polarized [I6]. At the metal/insulator interface, the screening of the 
polarization charge alters the FE polarization orientation resulting in a linear change 
of the surface magnetization. The resulting magnetization in the FM structure decays 
exponentially away from the interface. Taking the FM material as an ideal metal with 
the screening length of around [17] lA, the exchange interaction between the additional 
surface magnetization and the FM part is thus limited to only the first site. Switching 
to dimensionless units, we introduce the reduced polarization (pj(t) = Pj{t)/Ps) and 
magnetic moment (Sj(t) = /Xj(t)//is) vectors, where Ps is the spontaneous polarization 
of (bulk) BaTiOa and /xg is the magnetic moment at saturation of (bulk) Fe. 

2. Theoretical formalism 

The total energy of the FE/FM-system in a very general one- dimensional case consists 
of three parts 

Fj: = -Ffe + -^FM + Ec- (1) 
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The ferroelectric energy contribution reads 

^FE -a Fs 2_^ y-^ — Vj + -^- Pj + 

i=o 

^(p..i-p.f-PrE(t)), (2) 

whereas the ferromagnetic energy part is 

Ffm= 5^ (-J.S,-S,+i-A(5n'-/^sSi-B(t)j. (3) 

i=0 

E(t) and B(t) are respectively external electric and magnetic fields. 

Various pinning effects that may emerge in both the FE and the FM parts due to 

imperfections and expressed by kfej and Jj, Di respectively, are not considered here, 

i.e. KpEj = ^fe, Ji = J, Di = D. 

We consider that the linear FE/FM coupling appears as a result of the exchange 

interaction of the magnetization induced by the screening charge at the interface and the 

local magnetization in the ferromagnet (for details we refer to [18]) and can be written 

as 

Ec = APs/isPo-So. (4) 

The meaning of the quantities appearing in these equations is explained in Table [TJ 
Based on the parameters obtained from ab-initio calculations for BaTiOs/Fe- 
interface[T3l [T9]. we estimate the coupling constant as A = J Ofm'^s/Isfe £0/^0/^3) ^ 
2 ■ 10~^ s/F, where the surface ME-coupling constant is as = 2 ■ 10~^° G cm^/V. We find 
this value is too low to obtain a sizable ME-response for the one-dimensional multiferroic 
interface. In what follows, we vary A and explore the dependence of the multiferroic 
dynamics on it. 

The polarization dynamics is governed by the Landau-Khalatnikov (LKh) equationj20[ 
|2l],i.e. 

luPs-jf = Hj = — ^-—^ — , (5) 

where 'ju is the viscosity constant (Table [1]) and H • stands for the total external and 
internal fields acting on the local polarization. The magnetization dynamics obeys the 
Landau-Lifshitz-Gilbert [22] (LLG) equation of motion 

^ = -T-^ [S. X H™(t)] - -p^ [S. X [S. X H™(t)]] , (6) 
dt 1 + a'^FM 1 + a FM 

where 7 is a gyromagnetic ratio (Table [T]) and aFM is the Gilbert damping 

parameter. The total effective field acting on Sj is defined as a sum of 

1 f)I^ 

deterministic and stochastic parts Hf^(t) = — — h CA^)- The characteristics of 

yUs obi 
the additive white noise associated with the thermal energy k^T are [23] (dkit)) = 

and {Cik(t)Crni(t + At)) = 6imSkiS(At). Here i and m index the corresponding 

/^s7 
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sites in the FM-material. k and / are the Cartesian components of ( and At is the time 
interval. The coupled equations of motion ([5]) and ([6]) are solved numerically in reduced 
units, re-normalizing the energy ([T]) over doubled anisotropy strength D. Thus, the 
dimensionless time in both equations is r = coa t = 7 -Ba t = 7 2D t/^s and the reduced 
effective fields are h™(r) = B.™{t)/Ba, hf = Hf /(77,Ps5a). 
To endorse our results we conducted furthermore kinetic Monte Carlo (MC) simulations 
for the description of the dynamics of a FE/FM chain subjected to external magnetic 
and electric fields. The electric dipoles and the magnetic moments are understood as 
three-dimensional classical unit vectors, which are randomly updated with the standard 
Metropolis algorithm |24]. The period of the external field is chosen to be 600 MC steps 
per site. 
In the following the system is described via the reduced total polarization ps(^) 

AfpE-l 

j=0 

and the reduced net magnetization Ss(t) 

AfpM-l 

S^W = 1^ E S.W- (8) 

3. Results of numerical simulations 

To demonstrate the response of the FE/FM-chain to external fields using the LKh and 
the LLC equations we used a damping parameter apM = 0.5, which is significantly higher 
than the experimental value[25j (a(Fe)=0.002), in order to achieve a faster relaxation 
of the net magnetization for both B- and E-drivings. It was assured in our calculations 
that the FE-subsystem is far from its phase transition temperature at T = 10 K and the 
FM-subsystem remains non- superparamagnetic on the whole time scale of consideration. 
Fig. [2] shows the hysteresis loops in the presence of a harmonic external magnetic field 
i?z(t) = Bo^cosut. The amplitude of the field Bq^ is chosen to be comparable with 
the exchange interaction energy J as well as the coupling energy Ec- The period 
of the external magnetic field is chosen to exceed the field-free precessional period 
{27r/u ~ 52^pi"ec'j q£ ^]-^g LLC. Irrespective of the temperature and the calculational 
method this field is capable of switching the magnetization of the FM-chain (Fig. [2b, d). 
The FE-polarization indirectly driven by the external magnetic field is not completely 
switched according to both methods (Fig. [2^, c). The role of thermal fluctuations on 
the FM part only (cf. eq. ([5])) is exposed by the p(B)-behavior shown in Fig. [2^, 
whereas the MC method accounts for temperature effects also on the polarization (Fig. 
[2t). As a result, the p(i?)-hysteresis (Fig. [2t) shows a clear temperature dependence, 
which becomes especially pronounced for the one-dimensional chain in which thermal 
fluctuations degrade the polarization/magnetization ordering more intensively than for 
the case of a two-dimensional system. 
The hysteresis loops for the external electric field of the form Ez{t) = Eo^cosut are 
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presented in Fig. El The energy of the apphed electric field is comparable with the 
coupling energy. As a result, the total polarization can be completely switched (Fig. 
[3^, c). As inferred from Fig. |3b, d the net magnetization is not fully switched. Only 
several first spin sites follow the electric field due to the coupling at the interface. 

4. Discussion 

In real experiments several additional effects may affect the polarization and the 
magnetization dynamics. Below we estimate how strong these effects might be for 
the multiferroic chain. 

4-1- Effect of depolarizing fields 

In the general case the field E(t) entering equation ([2]) is an effective field that consists 
of the applied electric field, e.g. Ez{t)e^, and the internal depolarizing field Eqf created 
by the screening charges (SC) at the interface. 

Generally, the depolarizing field may well be sizable and affects thus the dynamics 
[301 ISTl [32] . Here we estimate the strength of Edf by introducing a one- dimensional 
SC at the interface Qsc = -Psct^ (Table HI). As a result, the electric field induced 
by the SC is opposite to the local polarization (Fig. [1]) and can be written as 
Edf = — Qsc/(4vr£:o£^FEa^^^)ez, where Eq = 8.85 • 10"^^ A-s/(V-m) is the permittivity of 
free space, 5fe ~ 2000 is the dielectric constant in barium titanate and Uj is the index 
numbering the polarization sites starting from the interface, e.g. nj=Q = 1. Thus, the 
strength of the depolarizing field calculated for Uj^q = 1 and upon the other parameters 
is Edp ^ 2 ■ 10^ V/m, which is at least one order of magnitude smaller than the 
amplitudes of the applied electric field (~ 4 ■ 10^ V/m). Keeping in mind that i?DF 
decays in the FE away from the interface we can neglect the depolarizing field. 

4-2. Effect of induced electric and magnetic fields 

According to the Maxwell's equations an oscillating magnetic field induces an oscillating 
electric field, and an alternating voltage produces an oscillating magnetic field. 
The situation becomes especially important for the second case, since even small 
induced magnetic fields aligned perpendicularly to the inducing field and hence to the 
initial state of the magnetization can sufficiently assist the switching at appropriate 
frequencies [33] l3lj . 

From the Faraday's law of the Maxwell's equations V x E = —/J-FEl^o^f and for the 
given applied magnetic field B^it) = Bq^ cos cot the induced electric field acting on the 
FE-polarization is oriented perpendicularly to the inducing field {XY-plane, Fig. [1]). 
Its amplitude for relative magnetic permittivity in BaTiOs [35] /xfe ~ 1 and /xfm ~ 5000 
scales as E^ff''^ = aFE/^FE/(/^FM) Bq^cj f» 2 V/m. 

Likewise, when an external electric field -Ez(t) = EozCosut is applied, according to the 
Ampere's law of the Maxwell's equations V x B = /xfm A^o^fm^o ^; the direction of 
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the induced field is perpendicular to the inducing field. The amplitude of the induced 
magnetic field in iron for ofm = 0.28 ■ 10~^m and epM ~ 1 (since iron is assumed as an 
ideal metal) is S"'^ = ofm A^fm /^o ^fm £o -E'oz w ~ 2 ■ 10~^ T. 

Our calculations with the estimated (and even larger) amplitudes of the induced electric 
field show no infiuence on the Z-projection of the FE-polarization. This is a consequence 
of the uncoupled nature for the projections of the FE-polarization (cf. equation ([5])). 
Both numerical methods give a slightly enhanced (less than 1%) ME-response in the 
presence of the induced magnetic field. Therefore, the effect of the induced electric and 
magnetic fields can be deemed irrelevant for the considered multiferroic chain and for 
the chosen range of frequencies. 

4.3. Frequency dependence of the magnetoelectric response 

A variation of the frequency u of the external electric and magnetic fields can also affect 
the ME-response of the multiferroic chain. 

Fig. m demonstrates the response of the multiferroic interface to an external magnetic 
field. The periods of the external magnetic field should be compared with a characteristic 
field-free precessional time Tp^^ ?a 4 ps which is valid for bulk iron. As one expects, 
magnetic fields with longer periods favor better saturation of the magnetization (Fig. 
Hb). The response of the total electric polarization to the external magnetic field 
becomes enhanced with increasing period of B-field. This is confirmed by both numerical 
methods (Fig. Hti, c). 

The multiferroic response to an external electric field is shown in Fig. [5] for which 
the situation of very short electric fields (less or around Tp^^) is addressed. The 
magnetization does not relax quick enough resulting in the form of a hysteresis which 
is similar to the ferroelectric one. Additionally, we obtain an increase of the net 
magnetization response (Fig. [5b). This feature can also be observed using the MC- 
method (Fig. [Sji). 

5. Summary 

The main result obtained using two independent methods - the direct solution of 
the LKh and the LLG equations ([51 [6]) and the kinetic MC method - is that due to 
the coupling at the interface of FE/FM the ferromagnetic subsystem responds to an 
external electric field and the ferroelectric subsystem responds to an external magnetic 
field. A use of both methods allowed a comparison of dynamical and statistical 
approaches for studying coupling phenomena at the FE/FM interface. Additionally, 
the LKh/LLG equations provide an insight into the real time temporal behavior, while 
MC approach is very useful to inspect the temperature infiuence on both sides of the 
interface. 

This research is supported by the research projects DFG SFB762 (Germany) and 
FONCICYT 94682 (Mexico). 
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FE (BaTi03) ]y ' FM (Fe) 



Pl Po So Si S2 S49 

Figure 1. Alignment of electric dipoles p, and magnetic moments S; in the considered 
one-dimensional chain. The directions of Si{t = 0) ~ {0, 0, 1} and p (i = 0) = {0, 0, 1} 
represent the initial configuration. The easy axis in the FM part is along the Z- 
direction. 



Table 1. Parameters used in the numerical calculations. 



FE-material (BaTiOs) 


Number of sites iVpE 




50 


Polarization 26] Ps, C/m^ 




0.265 


Initial state p^.(t=0), [Ps 




{0.0,0.0,1.0} 


Constant [20] 7^, [Vms/C 




2.5-10~5 


Constant [26] apE, [Vm/C 




-2.77-10^ 


Constant [26J /3fe, [VmVC^; 




1.70-10^ 


FE-interaction kfe, [Vm/C 




1.0-10^ 


Lattice constant j27j a. m 




0.4- 10-9 


FE/FM-coupling A, [s/F 




parameter 


FM-material 


(Fe) 




Number of sites A^fm 




50 


Gyromagn. ratio 7, [(Ts)""*^ 




1.76 ■ 10" 


Moment per site 28] /is- /^p 




2.2 


Initial state Sj(t=0), [/is 


{0 


14,0.14,0.98} 


Anisotropv strength j29 D. J 




1.0-10-22 


Exchange strength j29 J, J 




1.33-10-21 


Damping apu 




parameter 
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Figure 2. The reduced total polarization/magnetization response to external 
magnetic field of the form Bz{t) = Bq^ cos cot. The loops a) and b) are obtained 
by using the LKh and the LLG equations; c) and d) are calculated using the MC 
method. In the both methods parameters are chosen such that Ec ~ HsBz ~ J, i.e.: 
A = 240 s/F, Boz = 6.65Ba, oj = 3.61 • 10" s-\ Eq^ = V/m, a = 0.5. 20 first 
periods (1/w) are omitted; the hysteresis curves are averaged over 100 (a, b) and 200 
(c, d) subsequent periods. 
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Figure 3. Hysteresis loops of tlie total reduced polarization/magnetization as a 
function of external electric field of the form Ez{t) — Eo^cosojt. The curves a) and 
b) are obtained by using the LKh and the LLG equations; c) and d) are calculated 
using the MC method. Parameters are chosen such that Ec ~ a^PsE^ ~ J, i.e.: 
A = 240 s/F, Boz = T, Eq^ = 4.07 • 10^ V/m, uj = 3.61 ■ 10" s-\ a = 0.5. 20 first 
periods (1/w) are omitted; the hysteresis loops are averaged over 100 (a, b) and 200 
(c, d) subsequent periods. 
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Figure 4. Response of the one-dimensional multiferroic structure to the time- 
dependent magnetic field Bz{t) = Bo^coaivt plotted for various frequencies w. The 
field-free precessional time in the FM part is T^^'^ = T^'^^'^ « 4 ps. The curves a) and 
b) are obtained by using the LKh and the LLG equations; c) and d) are calculated using 
the MC method. Parameters are Tq = K, a == 0.5, Bqz = 6.65Ba and A = 240 s/F. 
Several first periods are omitted. 
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Figure 5. Response of the one-dimensional multiferroic structure to the time- 
dependent electric field £'z(i) = Eqz cos cot plotted for various frequencies w. The 
field-free precessional time in the FM part is Tp^'^ = T^'^'^'^ « 4 ps. The curves a) and 
b) are obtained by using the LKh and the LLG equations; c) and d) are calculated 
using the MC method. Parameters are Tq — K, a — 0.5, -Eqz — 4.07 • 10^ V/m and 
A — 240 s/F. Several first periods are omitted. 



